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CANGAROO-II Imaging Atmospheric Cherenkov Telescope (lACT). The flux 
was 0.12 times of that of Crab at 1 TeV. The signal centroid is consistent with 
the peak of the X-ray emission in the north-west rim of the remnant. 

Subject headings: gamma rays: observation — (stars:) supernovae: individual 
(RX J0852.0-4622) 

1. Introduction 

Young supernova remnants (SNRs) have for many years been believed to accelerate cos- 
mic rays up to energies around the "knee" of the cosmic ray spectrum at ~100TeV, however 
direct observational evidence has only recently been found to support this. The interactions 
of such high energy cosmic rays with surrounding material produce a signature in TeV en- 
ergy gamma-rays, and so TeV observations of SNR are therefore of fundamental importance 
in understanding the origin of Galactic cosmic rays. To date, the only SNRs which have 
been reported to emit TeV gamma-rays are the shell- type SNRs SN1006 (Tanimori et al. 
1998), RX J1713. 7-3946 (Muraishi et al. 2000; Enomoto et al. 2002b), and Cassiopeia A 
(Aharonian ct al. 2001). However, although RX J1713.7— 3946 has been detected by two 
groups (Hofmann 2004), there have as yet been no confirming reports for Cassiopeia A, and 
for SN1006, a recent upper limit is significantly below the originally reported fiux (Hofmann 
2004). Therefore, in addition to further observations of these sources, it is necessary to 
search for further examples of such SNRs. 

RX J0852.0-4622 (G266.2-1.2) is a SNR located (in projection) at the southeast corner 
of the Vela SNR. It was discovered by Aschenbach (1998) using data from the ROSAT All-Sky 
Survey. It is one of the few known SNRs which display strong non-thermal X-ray emission 
(Tsunemi et al. 2000; Slane et al. 2001). In addition, the 1157 keV fine of ^''Ti was detected 
by COMPTEL (lyudin et al. 1998). Prom a consideration of the hfetime of ^^Ti (~90yr), the 
X-ray flux, and the angular size of the X-ray emission, ~2°, the age and distance from the 
Earth were estimated to be 500-1100 yr and 80-500 pc, respectively (Aschenbach, lyudin, & 
Schonfcldcr 1999). An independent estimate of the distance, of 1-2 kpc, was provided by 
Slanc ct al. (2001) from a simple scaling of the column density, which suggested that the SNR 
was at least several times more distant than the Vela SNR. ASCA hard X-ray images revealed 
a shell-like morphology with a featureless spectrum well described by a power-law (Tsunemi 
et al. 2000; Slane et al. 2001). Radio emission also has a non-thermal spectrum (Combi, 
Romero, & Benagha 1999; Duncan & Green 2000). Both the X-ray and radio emission 
are brightest in the north-west rim of the SNR. These data suggest that particles may be 
being accelerated to ~100TeV. At these energies electrons can give rise to TeV gamma-rays 
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through inverse Compton scattering, and protons can interact with surrounding material to 
produce TeV gamma-rays via 7r° decay. 

Furthermore, RX J0852.0-4622 is, hke RX J1713.7-3946, an extended SNR at X-ray 
energies. The relatively large angular size enables us to study the detailed morphology even 
with the existing angular resolution. 

We have observed RX J0852.0-4622 with the 10 m diameter CANGAROO-II Imag- 
ing Atmospheric Cherenkov Telescope (lACT). The telescope (Kawachi et al. 2001) detects 
gamma-rays above several hundred GeV by reconstructing optical Cherenkov images gener- 
ated by relativistic secondary particles in the cascades produced when high-energy gamma- 
rays (and background cosmic rays) interact with the Earth's upper atmosphere. Parametriz- 
ing the optical images enables the direction and energy of gamma-ray events to be inferred. 
The telescope is located near Woomera, South Australia (136°47'E, 31°06'S), from which 
RX J0852. 0—4622 culminates at ~15° from the zenith, enabling observations with a lower 
energy threshold. The angular resolution of CANGAROO-II was estimated to be 0.30-0.39° 
(2.6-3.4 pc at a distance of 0.5 kpc) with an energy threshold of 500 GeV depending on the 
observing conditions and the spectrum. 



2. Observations 

The observations were carried out between 2001 December 12 and 2002 February 15 (19 
nights), and between 2003 January 5 and February 28 (25 nights). The telescope tracked the 
peak of the X-ray emission, in the north-west rim of the SNR {a — 132.25°, 6 — —45.65°, 
J2000 coordinates). The field of view (FOV) of the camera is 2.76° x 2. 76°. The brightest 
star in the FOV, SAO 220422, has a visual magnitude of 4.1, and this region is brighter 
than our typical on-source FOV. As a result, we used a higher trigger threshold, requiring 
5 triggered pixels rather than the usual 3 (Itoh et al. 2003). The recently installed lights 
several kilometers from the telescope were expected to have little effect on these observations 
as the telescope was pointed close to the zenith (Itoh et al. 2003), and this was confirmed 
from the azimuthal angle dependence of the shower rate (see §3). Each night was divided 
into two or three periods, i.e., ON-OFF, OFF-ON-OFF, or OFF-ON observations. ON- 
source observations were timed to contain the meridian passage of the target, as was done 
by Enomoto et al. (2002b). In total, 5900 min. of ON- and 5300 min. of OFF-source data 
were obtained. 
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3. Analysis 

First, 'cleaning' cuts on camera images were applied, requiring (0.115°-square) pixel 
pulse-heights of greater than 3.3 photoelectrons, and Cherenkov photon arrival times within 
±50 ns of the median arrival time. Clusters of at least five adjacent triggered pixels (rather 
than the usual four-pixel cut) were required in each event to minimize the effects of the 
bright star field. After these pre-selection cuts, which reduced events due to background 
light were reduced by 99%, the shower rate was stable on a run-to- run basis for observations 
in the same year. The systematic difference of the run-by-run acceptance within the same 
year is expected to be less than 12%. The ON/OFF shower rate differences in 2002 and 2003 
were —8 ± 15% and — 1 ± 13%, respectively. By examining the event rates within each run 
we were able to reject periods affected by cloud, dew forming on the mirrors, instrumental 
abnormalities, etc. Only data taken at elevation angles greater than 60° were accepted. 
After these cuts, 4300 min. of ON- and 3900 min. of OFF-source data survived. 

Trigger rates for each pixel per 700 fis were monitored by a scaler circuit in real-time 
and recorded each second. These data were used to exclude 'hot' pixels (generally due to 
the passage of a star through the FOV of a pixel) in off-line analysis. Hillas parameters were 
then calculated to discriminate gamma-rays from cosmic rays based on the image shape and 
orientation (Hillas 1985). Further, in each year's data, we masked a small number of pixels 
which showed deformed ADC spectra, possibly due to a hardware fault. Discrimination of 
the cosmic ray background from gamma-rays was carried out using the likelihood method of 
Enomoto et al. (2002a). 

4. Results 

The resulting distributions of the image orientation angle, a, for the combined data 
in 2002 and 2003 are shown in Fig. 1. The normalizations between the ON- and OFF- 
distributions were carried out using data with a > 27°. The numbers of excess events 
(a < 18°) were 530 ± 120 (in an observation time of 2100 min.), 540 ± 140 (2200 min.), 
and 1080 ± 180 (4300 min.) in 2002, 2003, and the combined data, respectively. The excess 
rates for 2002 and 2003 were similar to each other. Nightly signal rates were also checked 
during both years. The largest deviations occurred with a rate 2.5 ± 0.9 times larger than 
the average for 2002 and 3.0 ± 1.0 times for 2003, respectively. These are not unexpected 
statistically, and therefore there is no evidence of time-variability in the TeV emission. 

To check on the spatial distribution of the signal we derived the "significance map", 
shown by the blue contours in Fig. 2. The contours were calculated from the distribution of 
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the detection significance determined at each location from the difference in the a plots (ON- 
minus OFF-source histogram) divided by the statistical errors. The centroid is consistent 
with position of the X-ray maximum, within our possible systematic uncertainty of 0.1°. 
The region enclosed by the 65% (of the peak TeV significance) contour is elliptical with 
a semi-minor axis of 0.19° (NE-SW) and a semi-major axis of 0.35° (NW-SE), while our 
angular resolution was estimated to be 0.30-0.39°. The acceptance of the CANGAROO-II 
telescope decreases smoothly with offset from the tracking center, faUing to 50% at a 0.9° 
offset, i.e. at the edge of the trigger region. It is difficult to obtain the reliable upper limits 
on the two neighboring X-ray enhancements in the rim since they are near the edge of the 
trigger region. 

After correcting for this acceptance, the differential fluxes listed in Table 1 were derived. 
As both statistical and systematic errors are included, the energy bins overlap somewhat, par- 
ticularly at low energies. The systematic uncertainty for the energy determination (~20%) 
dominates the errors in the energies. The Spectral Energy Distributions (SEDs), derived 
from radio, X-ray, MeV gamma-ray, and our TeV gamma-ray observations, are plotted in 
Fig. 3. The hatched area indicates the allowed region, where the foi' the power-law fitting 
< x^jj^ + 1 including energy uncertainty, and x^^j^ is the for the best fit. The values are 
listed in Table 1. The derived spectrum has a power-law index of — 4.3l4;4, the errors of 
which come predominantly from the uncertainty in energy. 

The derived spectrum from RX J0852. 0—4622 seems rather steep. The photon index of 
—2.6, however, is within 1 a level due to the large uncertainty on the energy determination. 
This indicates that we cannot firmly conclude that the energy spectrum as soft as spectral 
index of —4.3. Recently CANGAROO-II collaboration has reported a similar steep spectrum 
from the Galactic Center (Tsuchiya et al. 2004). This spectrum is different from that reported 
by H.E.S.S. collaboration (Aharonian et al. 2004a). Applying the same procedure to the 
Galactic Center result, the allowed spectral index range became — 4.615;q. further answer 
on this is found in this stage. Time variability might be one of solution (Aharonian & 
Neronov 2004). On the other hand, the difference in the energy spectrum from the SNR 
RX J1713.7-3946 between the CANGAROO-11 and the H.E.S.S. is 2.2a level (Aharonian 
et al. 2004b). These complex situations encourage the further observations and analysis of 
the stereoscopic data using lACTs. 

Our Monte Carlo simulations predict, even for a point-source, that the a distribution of 
gamma-ray events for this soft spectrum should have a broader peak than that for a Crab-like 
spectrum. Therefore we cannot conclude from the experimental data whether this source is 
extended or point-like. 

Various checks on the signal yield and position were carried out, by varying thresh- 
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olds, clustering cuts, Hillas parameter values, etc: these yielded consistent fluxes within the 
systematic errors given in Table 1. Also, Crab nebula data were analyzed with the same 
code, with the derived flux and morphology consistent with previous measurements and the 
point-spread function, respectively. The distributions of the Hillas parameters for the ex- 
cess events were checked and found to be to be consistent with Monte Carlo simulations for 
gamma-rays. The OFF-source data were compared with the Monte Carlo simulations for 
protons, and found to be consistent. 



5. Discussion 

The gamma-rays around 1 TeV are most likely produced by one of two mechanisms. We 
consider flrst the synchrotron/inverse Compton (IC) model. The spectra obtained using this 
model are shown in Fig. 3 (upper). This model has the following parameters: the magnetic 
field {B), the spectral index of electrons (7), and the maximum electron energy (-Ee,max)- The 
cutoff energy of the synchrotron emission is proportional to E^^^^B, which was constrained 
by the X-ray and the radio data. The cutoff of the IC emission is proportional to -Be.max- 
The ratio of the peak of the synchrotron component and the peak of the IC component in 
the SED is proportional to B^. li B is too small, -Eg max is too high to explain our data, as 
shown by (a) in Fig 3. The lower limit of B in this model is ~12 /xG as shown by (b) in 
Fig. 3. This model cannot explain our data. 

A more complex model in which the synchrotron and IC emissions come from different 
zones (a two-zone model) was considered by Aharonian, Atoyan, & Kifune (1997). At first, 
we assume the same electron spectra and B with different zone sizes. The lines(c) and (d) 
in Fig. 3 show the IC emissions with size ratios (Vxcv/Vx-ray) of 1 and ~ 10^, respectively. 
The ffux ratio between lines (c) and (d) corresponds roughly to the size ratio. This model 
requires a high size ratio and a strong magnetic field (~ 1.6 mG). Chandra observations 
revealed small-scale structures in the rim of SN1006 (Bamba et al. 2003): similar high 
resolution observations (in both X-rays and TeV gamma-rays) of RX J0852. 0—4622 will 
help constrain the size ratio and refine this two-zone model. Stereo observations in the TeV 
region are required. 

The other possible origin for the gamma-rays is 7r° decay. The estimated fiux from 
7r° decays is shown by the red line in Fig. 3 (lower). The emissions from the electrons 
were also considered using the electron/proton ratio (Kcp). Bremsstrahlung emission, shown 
by the green line, is constrained by the ASCA data. Therefore Kep was limited to be at 
most 4x10^''. As a result, both the ffux of the inverse Compton emission and that of the 
bremsstrahlung were too small at TeV energies to explain our data. The sum of the tt^ 
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emission, the inverse Compton emission and the bremsstrahlung is shown by the hght blue 
hne. The spectrum of accelerated protons is assumed to be a power-law with an exponential 
cutoff. The 7r° model has the following parameters: the spectral index of the protons, the 
maximum energy of protons (-E'p,max)) and A = (£'o/10^°6rg)(no/protons cm^^)((i/0.5kpc)^^, 
respectively. Here, Eo, no, and d are the total energy of protons in the region from which 
gamma-rays were detected, the number density of protons where the interactions occur, and 
the distance from the earth, respectively. The spectral index and the maximum energy of 
protons were assumed to be those of electrons, i.e., —2.5 and 8TeV, respectively. For the 
best-fit spectrum A =4.9 was obtained. Assuming the distribution of the energy of the 
accelerated cosmic rays is isotropic, the total energy of accelerated cosmic rays in the whole 
SNR is proportional to the ratio Vtot/VTeV, where Vtot is the total volume. This ratio was 
estimated to be 46 from the "significance map" with the assumption of spherical symmetry. 
The range of -Eo(Vtot/VTcv) of 10^^-10^'^ergs corresponds to Uq of 23000-230, which are 
typical values for molecular clouds. Thus, the model of 7r° decays can naturally explain the 
multi-band spectrum. 

The molecular distribution toward the whole extent of the Vela SNR has been mapped 
by Moriguchi et al. (2001). As RX J0852. 0—4622 is embedded (in projection) in the Vela 
SNR, it is not possible to uniquely associate the detected molecular clouds with one or other 
SNR. The general anti-correlation of the X-rays and ^^CO with the Vela SNR suggests a 
strong interaction between this SNR and the interstellar medium, but does not rule out the 
existence of molecular clouds around RX J0852. 0—4622, encouraging the further analysis. 
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a(degree) 

Fig. 1. — Distributions of a (image orientation angle) for the combined data in 2002 and 
2003. The points with error bars show the ON-source data and the hatched histogram is 
the OFF-source data, which was normahzed to the ON-source data using the number of the 
events with a > 27°. 
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Fig. 2. — The "significance map" obtained by the CANGAROO-II telescope is shown by the 
blue contours. The telescope tracking center is at the north-west rim of RX J0852. 0—4622. 
The red contours with levels at 20%, 45%, 65%, and 80% are an ASCA GIS image. The 
green contours show the 4850 MHz radio emission. 
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Fig. 3. — Spectral energy distribution of the north-west rim of RX J0852. 0—4622. The 
points with error bars and the arrows in the TeV region are the fluxes and the 2a upper 
limits, respectively, from our observations. The hatched area indicates the allowed region, 
where the for the power-law fitting < x^jj^ + 1 including energy uncertainty, and Xmin 
is the for the best fit. The arrows around 100 MeV show the la upper limits obtained 
using the EGRET data of the diffuse emission. They were derived from the count rate plus 
twice the standard deviation based on the archival data of the EGRET gamma-ray intensity 
map, which was integrated over the area considering the angular resolutions of EGRET and 
CANGAROO. The open squares with the error bars show the Parkes radio data (Duncan 
& Green 2000). The black filled areas show the ASCA X-ray data (Slane et al. 2001). The 
red line shows the estimation for vr^ gamma-rays, the details of which are given in the text. 
The blue lines shows the spectra of the various synchrotron/inverse Compton (IC) models. 
The green line shows the bremsstrahlung emission from electrons. The light blue line shows 
the sum of the vr", the inverse Compton, and the bremsstrahlung emissions. 
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Table 1: Differential fluxes. 



Mean energy Flux 

of bin [TeV] [ph/cmVs/TeV] 



Lower bound* Upper bound* 



0.48±0.11 
0.55±0.12 
0.79±0.18 
1.93±0.48 

5.84±1.51 
10.38±3.15 



(9.43±3.49)xl0-^^ 

(3.85±1.78)xl0-^^ 
(8.54±4.63)x 10-^2 

(2.80±2.45)x 10-^3 

3.18x10"^^ ** 
1.91x10-14 ** 



1.99x10"" 
1.03x10-11 
1.47x10-12 

1.35x10-1^ 



3.01x10-1° 
1.28x10-1° 
1.84x10-11 
1.22x10-12 



* These were estimated where the for the power-law fitting < xLin + 1 including energy uncertainty, where 

Xmin ^^'^ ^^St fit. 

** These are 2 a upper limits. 



